
Claims : 

1 1. A method of processing a digital image, comprising the 

2 steps of: 

3 providing digital data indexed to represent positions on 

4 a display, said digital data being indicative of an intensity 

5 value I-(x,y) for each position (x,y) in each i-th spectral 

6 band ; 

7 adjusting said intensity value for said each position in 

8 each i-th spectral band to generate an adjusted intensity 

9 value for said each position in each i-th spectral band in 
W accordance with 

§ J2^ n dog JjOcy) - log [l d (x,y) *F n (x,y)]) , i = l,...,S 

fe where S is the number of unique spectral bands included in 

P 

12 said digital data and, for each n, W n is a weighting factor 

S 

1TB and F n (x,y) is a unique surround function applied to said each 

ill 

fi£ position (x,y) and N is the total number of unique surround 

6 

$5 functions; 

16 filtering said adjusted intensity value for said each 

17 position of said image in each of said S spectral bands using 

18 a filter function wherein a filtered intensity value R f (x,y) 

19 is defined; and 

20 selecting a maximum intensity value V^x^) from the 

21 group consisting of said intensity value I i (x,y) and said 

22 filtered intensity value R^x^). 
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2. A method according to claim 1 wherein each said unique 
surround function is a Gaussian function. 

3, A method according to claim 2 wherein said Gaussian 
function is of the form 

-r 2 

satisfying the relationship 



unique constant for each of said N unique surround functions. 




where 




and, for each n, k n is a normalization constant and c n is a 
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4. A method according to claim 1 further comprising the step 
of multiplying said filtered intensity value R,- (x,y) by 



log 



Bl^x.y) 



1 = 1 



to define a color-restored intensity value R'^x^), where B 
is a constant, wherein said step of selecting using said 
filtered intensity value R-(x,y) is replaced with the step of 
selecting a maximum intensity value V^(x,y) from the group 
consisting of said intensity value I^x^y) and said color- 
restored intensity value R'^Xjy). 

5. A method according to claim 1 wherein said each position 
(x,y) defines a pixel of said display. 

6. A method according to claim 1 wherein, for each n, W n =l/N. 

CT: A method according to claim 1 further comprising the step 
of displaying an improved image using said maximum intensity 
value V 1 - (x, y) . 
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8. A method according to. claim 4 further comprising the step 
of displaying an improved image using said maximum intensity 
value V- (x, y) . 
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1 9 - A method according to claim 1 further comprising the steps 

2 of: 

3 defining a classification of said image based on dynamic 

4 range of said image in each of said S spectral bands; and 

5 selecting said filter function based on said 

6 classification of said image. 

1 10. A method according to claim 9 wherein said step of 

2 defining comprises the step of using image statistics 
Q associated with said image in each of said S spectral bands. 

is. 

%£ 

%&£ 
ft 

Si 11. A method according to claim 10 wherein said image 

(ja 

0 statistics include brightness and contrast of said image in 

a '3 each of said S spectral bands. 

ft 

ygL 12. A method according to claim 4 further comprising the 
steps of: 

3 defining a classification of said image based on dynamic 

4 range of said image in each of said S spectral bands; and 

5 selecting said filter function based on said 

6 classification of said image. 

1 13. A method according to claim 12 wherein said step of 

2 defining comprises the step of using image statistics 

3 associated with said image in each of said S spectral bands. 
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14. A method according to claim 13 wherein said image 
statistics include brightness and contrast of said image in 
each of said S spectral bands. 




15. A method of processing a digital image, comprising the 
steps of: 

providing digital data indexed to represent the positions 
of a plurality of pixels of a J-row by K-column display, said 
digital data being indicative of an intensity value I(x,y) for 
each of said plurality of pixels where x is an index of a 
position in the J-th row of said display and y is an index of 
a position in the K-th column of said display wherein a JxK 
image is defined; 

convolving said digital data associated with each of said 
plurality of pixels with a function 

-r 2 

a 2 

e c 

to form a discrete convolution value for each of said 
plurality of pixels, said function satisfying the relationship 



k is a normalization constant and c is a constant; 

converting, for each of said plurality of pixels, said 
discrete convolution value into the logarithm domain; 

converting, for each of said plurality of pixels, said 
intensity value into the logarithm domain; 



-i 




where 



r = y/x 2 +y 2 



subtracting, for each of said plurality of pixels, said 
discrete convolution value so-converted into the logarithm 
domain from said intensity value so-converted into the 
logarithm domain, wherein an adjusted intensity value is 
generated for each of said plurality of pixels; 

filtering said adjusted intensity value for each of said 
plurality of pixels with a filter function wherein a filtered 
intensity value R(x,y) is defined; and 

selecting, for each of said plurality of pixels, a 
maximum intensity value V(x,y) from the group consisting of 
said intensity value I(x,y) and said filtered intensity value 
R(x,y) . 

16. A method according to claim 15 wherein the value of said 
constant c is selected to be in the range of approximately 
0.01 to approximately 0,5 of the larger of J and K. 

17. A method according to claim 15 further comprising the 
step of displaying an improved image using said maximum 
intensity value V(x,y) . 

18. A method according to claim 15 wherein said step of 
filtering includes the step of selecting said filter function 
based on dynamic range of said JxK image. 
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19. A method of processing a digital image, comprising the 
steps of: 

providing digital data indexed to represent the positions 
of a plurality of pixels of an J-row by K-column display, said 
digital data being indicative of an intensity value I-(x,y) 
for each i-th spectral band of S spectral bands for each of 
said plurality of pixels where x is an index of a position in 
the J-th row of said display and y is an index of a position 
in the K-th column of said display wherein a (JxK),- image is 
defined for each of said S spectral bands and a JxK image is 
defined across all of said S spectral bands; 

convolving said digital data associated with each of said 
plurality of pixels in each i-th spectral band with a function 

-r 2 

e * 

for n=2 to N to form N convolution values for each of said 
plurality of pixels in each said i-th spectral band, said 
function satisfying the relationship 

-r 2 

k n J J e c?n dxdy = 1 

where 

r = \Jx 2 +y 2 

and, for each n, k n is a normalization constant and c n . is a 
unique constant; 

converting, for each of said plurality of pixels in each 



said i-th spectral band, each of said N convolution values 
into the logarithm domain; 

converting, for each of said plurality of pixels in each 
said i-th spectral band, said intensity value into the 
logarithm domain; 

subtracting, for each of said plurality of pixels in each 
said i-th spectral band, each of said N convolution values so- 
converted into the logarithm domain from said intensity value 
so-converted into the logarithm domain, wherein an adjusted 
intensity value is generated for each of said plurality of 
pixels in each said i-th spectral band based on each of said 
N convolution values; 

forming a weighted sum for each of said plurality of 
pixels in each said i-th spectral band using said adjusted 
intensity values; 

filtering said weighted sum for each of said plurality of 
pixels in each said i-th spectral band with a filter function 
wherein a filtered intensity value (x,y) is defined; and 

selecting a maximum intensity value (x,y) from the 
group consisting of said intensity value Ii(x,y) and said 
filtered intensity value R^(x,y). 

20. A method according to claim 19 wherein the value for each 
said unique constant c n is selected to be in the range of 
approximately 0.01 to approximately 0.5 of the larger of J and 
K. 
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21. A method according to claim 19 further comprising the 
step of multiplying said filtered intensity value R^x^y) by 



log 



Bl^x.y) 



2 = 1 



to define a color-restored intensity value R'^x^y), where B 
is a constant and S is a whole number greater than or equal to 
2, wherein said step of selecting using said filtered 
intensity value Rj (x,y) is replaced with the step of selecting 
a maximum intensity value V^(x,y) from the group consisting of 
said intensity value I^x^y) and said color-restored intensity 
value R' 1 - (x,y) . 

22. A method according to * claim 19 further comprising the 
step of displaying an improved image using said maximum 
intensity value V^x^y). 
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1 2 3. A method according to claim 21 further comprising the 

2 step of displaying an improved image using said maximum 

3 intensity value V-(x,y). 

1 24. A method according to claim 19 further comprising the 

2 steps of : 

3 defining a classification of said JxK image based on 

4 dynamic range of each said ( JxK) , image; and 

5 selecting said filter function based on said 

*6 classification of said JxK image. 

yQ 

% 2 5. A method according to claim 2 4 wherein said step of 
defining comprises the step of using image statistics 

?3 associated with each said (JxK)- image. 

0 

P 

[U 2 6. A method according to claim 2 5 wherein said image 

Ul 

£g statistics include brightness and contrast of each said ( JxK) 1 - 

U ' ■ ' 

3 image. 

1 27. A method according to claim 21 further comprising the 

2 steps of: 

3 defining a classification of said JxK image based on 

4 dynamic range of each said (JxK) ^ image; and 

5 selecting said filter function based on said 

6 classification of said JxK image. 
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28. A method according to claim 27 wherein said step of 
defining comprises the step of using image statistics 
associated with each said ( JxK) 1 - image, 

29. A method according to claim 28 wherein said image 
statistics include brightness and contrast of each said (JxK)- 
image . 



30. A method of processing a digital image, comprising the 
steps of: 

providing digital data indexed to represent positions of 
an image having S spectral bands for simultaneous output on a 
display, said digital data being indicative of an intensity 
value I s (x,y) for each position (x,y) in each i-th spectral 
band; 

adjusting said intensity value for said each position in 
each i-th spectral band to generate an adjusted intensity 
value for said each position in each i-th spectral band in 
accordance with 



where S is a whole number greater than or equal to 2 and 
defines the total number of spectral bands included in said 
digital data and, for each n, W n is a weighting factor and 



N 



Y, w n dog I i (x t y) - \oq[I ± (x,y) *F n (x,y) ] ) , i = 



F n (x,y) is a unique surround function of the form 



- r 



e 



satisfying the relationship 




where 



r = \Jx 2 +y 2 



and, for each n, k n is a normalization constant and c n is a 





unique constant where N is the total number of unique surround 
functions; 

filtering said adjusted intensity value for said each 
position in each i-th spectral band with a filter function 
wherein a filtered intensity value R^x^y) is defined; 

multiplying said filtered intensity value R-(x,y) by 



is a constant; and 

selecting a maximum intensity value V.(x f y) from the 
group consisting of said intensity value I^x.y) and said 
color-restored intensity value R^-Cx^y). 

31. A method according to claim 30 wherein, for each n, 
W n =l/N. 

32. A method according to claim 30 wherein the value for each 
said unique constant c n is selected to be in the range of 
approximately 0.01 to approximately 0.5 of the larger of J and 
K. 



log 



Bl^x.y) 



s 



£ Ji(x,y) 



to define a color-restored intensity value R^- (x,y), where B 



33. A method according to claim 30 further comprising the 
step of displaying an improved image using said maximum 
intensity value V,. (x,y). 



34. A method according to claim 30 further comprising the 
steps of: . 

defining a classification of said image based on dynamic 
range of said image in each of said S spectral bands; and 

selecting said filter function based on said 
classification of said image, 

35. A method according to claim 34 wherein said step of 
defining comprises the step of using image statistics 
associated with said image in each of said S spectral bands. 

36. A method according to claim 3 5 wherein said image 
statistics include brightness and contrast of said image in 
each of said S spectral bands. 



